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Abstract. We show that hadron production in relativistic heavy ion collisions at 
transverse momenta larger than 2 GeV/c can be explained by the competition of 
two different hadronization mechanisms. Above 5 GeV/c hadron production can 
be described by fragmentation of partons that are created perturbatively. Below 
5 GeV/c recombination of partons from the dense and hot fireball dominates. 
This can explain some of the surprising features of RHIC data like the constant 
baryon-to-meson ratio of about one and the small nuclear suppression for baryons 
between 2 to 4 GeV/c. 
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The Relativistic Heavy Ion collider (RHIC) has provided exciting data about 
hadron production at transverse momenta of a few GeV/c in central Au+Au collisions. 
The production of pions at high Pt was found to be suppressed compared to the scaled 
yield from p + p collisions 1 . This jet quenching effect can be understood by final 
state interaction of fast partons with the dense and hot medium produced in central 
heavy ion collisions. Fast partons lose energy via induced bremsstrahlung before they 
can fragment into high Pt hadrons . The suppression effect is dramatic and can be 
as large as a factor of 6 above Pt = 5 GeV/c. 

On the other hand the suppression of protons and antiprotons seems to be much 
less Experimental data from PHENIX show a proton/pion ratio of about 1 between 
1.5 GeV/c and 4 GeV/c 0]. This is surprising since the production of protons and 
antiprotons is usually suppressed compared to the production of pions because of the 
much larger mass. At high transverse momentum this can be understood in terms 
of perturbative quantum chromodynamics (pQCD) The fragmentation functions 
D a ^h{z) describe the probability for a parton a with momentum p to turn into a 
hadron with momentum zp, < z < 1. These fragmentation functions were measured 
for pions and protons, mainly in e + e~ collisions, and give a ratio p/tt° < 0.2 for 
Pt > 2 GeV/c when used in p + p and N + N collisions. The energy loss of partons 
in a medium can be taken into account by a rescaling of the parton momentum 6 . 
However this should affect baryons and mesons in the same way. 
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The lack of nuclear suppression for baryons is a challenge for our understanding 
of hadron production. The currently accepted picture assumes that a parton 
with large transverse momentum is produced in a hard scattering reaction between 
initial partons, propagates through the surrounding hot matter and loses energy by 
interactions, and finally hadronizes. Apparently the creation and interaction of a 
parton will happen independently of its later fate during hadronization. Therefore any 
unusual behavior that is different between baryons and mesons can only be attributed 
to hadronization itself. We propose to use recombination of quarks from the surface 
of the hot fireball as an alternative hadronization mechanism 0. 

In the fragmentation process a parton with transverse momentum pt is leaving the 
interaction zone while still being connected with other partons by a color string. The 
breaking of the string creates quark antiquark pairs which finally turn into hadrons. 
The distribution of one of these hadrons, which is bound to have less transverse 
momentum Pt = zpr, is then described by a fragmentation function. The average 
value (z) is about 0.5 for pions in p + p collisions. In other words the production of a, 
say, 5 GeV/c pion has to start with a 10 GeV/c parton in average, which are rare to 
find due to the steeply falling parton spectrum. Jet quenching even enhances the lack 
of high pt partons. On the other hand, the 5 GeV/c pion could be produced by the 
recombination of a quark and an antiquark with about 2.5 GeV/c each in average. 2.5 
GeV/c and 10 GeV/c are separated by orders of magnitude in the parton spectrum. 
The price to pay is of course that two of these partons have to be found close to each 
other in phase space. However we do have a densely populated phase space in central 
heavy ion collisions at RHIC where we even expect the existence of a thermalized 
quark gluon plasma. 

Recombination of quarks has been considered before in hadron collisions [5] and 
was also applied to heavy ion collisions §5 1. In QCD the leading particle effect in 
the forward region of a hadron collision is well known. This is the phenomenon that 
the production of hadrons that share valence partons with the beam hadrons are 
favored in forward direction. It has been realized that this can only be explained by 
recombination. This has fueled a series of theoretical work, see e.g. JTO] and references 
therein. Recently the recombination idea for heavy ion collisions, stimulated by the 
RHIC results, has been revived for elliptic flow and hadron spectra and ratios 

The formalism of recombination has already been developed in a covariant setup 
for the process of baryons coalescing into light nuclei and clusters in nuclear collisions 
jl.'il HI) . We give a brief derivation for the case of mesons. By introducing the density 
matrix p for the system of partons, the number of quark-antiquark states that will be 
measured as mesons is given by 

/J3 p 
(M;P\p ab \M;P) (1) 

Here \M; P) is a meson state with momentum P and the sum is over all combinations 
of quantum numbers — flavor, helicity and color — of valence partons that contribute 
to the given meson M. This can be cast in covariant form using a hypersurface £ for 
hadronization |15) 

„dN M _ fdaP-u(a) f d 3 q ( P \ w . . ( P \ 

s 
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Here E is the energy of the four vector P, da a measure on £ and u(a) is the future 
oriented unit vector orthogonal to the hypersurface £. w a and lOj, are the phase space 
densities for the two partons a and b, Cm is a degeneracy factor and 



*&(q)= y rfV*E(r,q) (3) 

where <I>]^ (r, q) is the Wigner function of the meson Jl] . 

Since the hadron structure is best known in a light cone frame, we write the 
integral over q in terms of light cone coordinates in a frame where the hadron has no 
transverse momentum but a large light cone component P + . This can be achieved 
by a simple rotation from the lab frame. Introducing the momentum k = P/2 — q of 
parton a in this frame we have k + = xP + with < x < 1. We make an ansatz for 
the Wigner function of the meson in terms of light cone wave functions <j>M (x) ■ The 
final result can be written as [7] 

l 

E iih = °v I da ^r^r I dxw a (a; X p + ) |<m*)| 2 (i - *)p + ). (4) 



d 3 P 1V1 J (2tt) 3 

s o 

For a baryon with valence partons a, b and c we obtain 

1 

E ~~jTj5 = CB J da ~T^W J dxidx 2 dxz5(xi+x 2 +x s ,-l) (5) 
s o 
x w a (a; x 1 P + )w b (a; x 2 P + )w c (a; x 3 P + ) \4>b(xi, x 2 , x 3 )\ 2 . 

<^s(xi, x 2 , X3) is the effective wave function of the baryon in light cone coordinates. 

A priori these wave functions are not equal to the light cone wave functions used 
in exclusive processes. We are recombining effective quarks in a thermal medium 
and not perturbative partons in an exclusive process. Nevertheless, as an ansatz for a 
realistic wave function one can adopt the asymptotic form of the light cone distribution 
amplitudes 

(j} M {x) = V30»(l - a), (6) 

<Pb(xi, £2, #3) = 12V35xia;2a;3 

as a model. However, it turns out that for a purely exponential spectrum the shape 
of the wave function does not matter. In that case the dependence on x drops out of 
the product of phase space densities 

w a (a;xP + )w b (a; (1 - x)P+) ~ e -^ + /^ e -(i-x)P+/T = e -p+/T (?) 

One can show that a narrow width approximation, using S peaked wave functions that 
distribute the momentum of the hadron equally among the valence quarks — 1/2 in 
the case of a meson and 1 /3 in the case of a baryon — deviates by less then 20% from 
a calculation using the wave functions in iffijl. That there is a small deviation can be 
attributed to the violation of energy conservation. Since recombination is a 2 — ► 1 
or 3 — > 1 process, energy will generally not be conserved if we enforce momentum 
conservation and a mass shell condition for all particles. However, one can show that 
violations of energy conservation are of the order of the effective quark masses or Aqce> 
and can therefore be neglected for transverse hadron momenta larger than 2 GeV/c. 
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Fragmentation of partons is given by 

° 

The sum runs over all parton species a and o~ a is the cross section for the production of 
parton a with momentum P a = P/z. We use a leading order (LO) pQCD calculation 
of a a J7| together with LO KKP fragmentation functions ^8]. Energy loss of the 
partons is taken into account by a shift of the parton spectrum by 

Apr = \pkpr- (9) 

We summarize that the transverse momentum dependent yield of mesons from 
recombination can be written as ~ Cmw 2 (Pt/2) m the simple narrow width 
approximation, whereas from fragmentation we expect ~ D{z) ® w(Pt/z). For an 
exponential parton spectrum w = e~ Pr / T the ratio of recombination to fragmentation 
is 

* = °»e-*(*-&) (10) 
F (D) y ' 

where (D) < 1 and (z) < 1 are average values of the fragmentation function and 
the scaling variable. Therefore R/F > 1. In fact, recombination always wins 
over fragmentation from an exponential spectrum (as long as the exponential is not 
suppressed by small fugacity factors). The same is true in the case of baryons. 

Now let us consider a power law spectrum w = A(Pt / fj)~ a with a scale fi and 
a > 0. Then the ratio of recombination over fragmentation is 

and fragmentation ultimately has to win at high Pt- We also note that we can expect 
a constant baryon/meson ratio from recombination, when the parton spectrum is 
exponential. In this case the ratio is only determined by the degeneracy factors 

dN B = Cb_ (l2) 
dN M C M ' [ 1 

For our numerical studies we consider fragmentation of perturbative partons and 
recombination from a thermal phase 

Wth = e -PTCOsHr,-y)/Te-y 2 / 2A2 ^3) 

with an effective temperature T. rj is the space-time rapidity and A«2 the width of 
the rapidity distribution. We fix the hadronization hypersurface S by the condition 
Tf = \lt 2 — z 2 — const. . We set 77 = 5 fm/c. A two phase parton spectrum with 
a perturbative power law tail and an exponential part at low transverse momentum 
is also predicted by parton cascades like VNI/BMS |T5] . 

The parameters, A for the average energy loss and T for the slope of the 
exponential parton spectrum, are determined by a fit to the inclusive charged hadron 
spectrum measured by PHENIX gOj. We obtain Awl GeV and T » 350 MeV. The 
temperature contains the effect of a blue-shift because of the strong radial flow. An 
additional normalization factor of about 1 /30 for the recombination part is necessary 
to describe the data. This is due to the use of an effective temperature which gives a 
too large particle number compared to the physical temperature that we expect to be 
around 175 MeV. 
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Figure 1. Left: the charged hadron spectrum in Au+Au collisions at yS = 
200 GeV as a function of Px- Contributions from fragmentation (dashed), 
recombination (dotted) and the sum of both (solid line) are shown. Data are 
from the PHENIX collaboration 1201 . Right: the spectrum of up quarks. 



In Figure H we show the inclusive charged hadron spectrum using recombination 
and fragmentation for pions, kaons, protons and antiprotons. We note that the 
crossover between the recombination dominated and fragmentation dominated regimes 
is around 5 GeV/c. It is in fact earlier for pions (~4 GeV/c) than for protons (~6 
GeV/c). In Figure^ we also give the spectrum of up quarks as an example for the 
partonic input of our calculation. We note that the crossover between the perturbative 
and the thermal domain is here around 3 to 3.5 GeV/c. It is characteristic for 
recombination that features of the parton spectrum are pushed to higher transverse 
momentum in the hadron spectrum. 

In Figure [21 we give the proton/pion ratio and the nuclear modification factor 
Raa for pions, protons and charged hadrons. The proton/pion ratio shows a plateau 
between 2 and 4 GeV/c in accordance with experiment and a steep decrease beyond 
that. This decrease has not yet been seen and is a prediction of our work. 

In the quantity Raa the effect of jet quenching is manifest for all hadrons at large 
transverse momentum. For pions Raa grows only moderately at low Pt in accordance 
with PHENIX data |21| . On the other hand recombination is much more important 
for protons, because they are suppressed in the fragmentation process. Therefore Raa 
nearly reaches a value of 1 below 4 GeV/c for protons. 

In summary, we have discussed recombination as a possible hadronization 
mechanism in heavy ion collisions. We have shown that recombination can dominate 
over fragmentation up to transverse momenta of 5 GeV /c. Recombination provides 
a natural explanation for the baryon/meson ratio and the nuclear suppression factors 
observed at RHIC. 
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Figure 2. Left: the p/w + ratio exhibits a transition between 4 and 6 GeV from 
the recombination dominated regime (p/n + ~ 1) to the fragmentation dominated 
regime (p/7r+ ~ 0.1). Right: the nuclear modification factor Raa f° r pions 
(solid), charged hadrons (dashed) and protons (dotted line). Data points are 
R AA for 7T° from PHENIX 
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